Measuring anisotropic scattering in the cuprates 
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A simple model of anisotropic scattering in a quasi two-dimensional metal is studied. Its simplicity 
allows an analytic calculation of transport properties using the Boltzmann equation and relaxation 
time approximation. We argue that the c-axis magnetoresistance provides the key test of this model 
of transport. We compare this model with experiments on overdoped Tl-2201 and find reasonable 
agreement using only weak scattering anisotropy. We argue that optimally doped Tl-2201 should 
show strong angular- dependent magnetoresistance within this model and would provide a robust 
way of determining the in-plane scattering anisotropy in the cuprates. 



I. INTRODUCTION 



The unusual transport properties in the normal state 
of the cuprates continue to attract widespread inter- 
est. Many aspects differ significantly from transport in 
conventional metalsu. The resistivity is linear in tem- 
perature) 2 ] with scattering which is apparently electronic 
rather than electron-phonon in originE The Hall effect 
shows strong temperature dependence^. On doping with 
zinc, both the resistivity and inverse Hall angle show a 
rather simple Matthiessen's rulecl, suggestive of additive 
scattering mechanisms. However, their temperature de- 
pendences differ significantly. This has led to the sug- 
gestion that the resistivity and the inverse Hall angle are 
controlled by very different scattering rates (the so-called 
'two- lifetime' behavior )n. This.-j/iew seems to be con- 
firmed by optical measurements^ which suggest two in- 
dependent scattering rates in these materials. The mag- 
netoresistance does not follow Kohler scaling behavior 
but appears also to be controlled by the Hall angle scat- 
tering rated. Attempts to understand this behavior may 
be divided into three categories. 

Anderson! has argued that a non-Fermi liquid de- 
scription must be invoked to understand transport in 
the cuprates. In Anderson's picture the electron decays 
into holons and spinons which separately control resistiv- 
ity and magneto-transport respectively. A phenomeno- 
logical transport equation was introduced by Coleman, 
Schofield and TsvelikS to capture a model of transport 
where the two lifetimes are controlled by independent 
fluids of particles. 

A second approach argues that the magnetic field plays 
a special role in the transport-process. In the model of 
Kotliar, Sengupta and Varmatl a singular skew scatter^ 
ing term is invoked. In the picture by Lee and LeetHI 
only recombined slave-fermion and slave-boson particles 
can interact with the true applied magnetic field. The 
slave particles themselves are insensitive to the true field 
because of the large fluctuations in the fictitious gauge 
field. In both of these scenarios the measured cyclotron 
frequency would appear to be temperature dependent — a 



result at odds with current optical Hall experimental]. 

The third possibility — the 'anisotropic scattering' 
scenario — envisages a conventional metallic state charac- 
terized by an electron quasiparticle scattering rate which 
depends strongly on the particle mnpentum. This pic- 
ture was first suggested by Cooped! 2 ] and has, been fur- 
ther elaborated upon by Stojkmdc and Pinesliia, and also 
by Yakovenko and coworkersLjEj. Pines et al. argue 
that antifcrromagnetic fluctuations in the normal state 
strongly scatter regions on the Fermi surface which are 
connected by the antiferromagnetic wave vector. Regions 
away from these 'hot-spots' are weakly scattered. How- 
ever, Hlubina and RiceES have claimed that the cold re- 
gions would tend to 'short-circuit' the hot spots, thereby 
leading to conventional transport at low temperatures. 
Ioffe and MillisEZl have introduced a new variant on the 
anisotropic scattering picture which is partly motivated 
by the short-circuiting problem and also inspired by the 
direct measurements of electron lifetimes. They use 
angle-resolved photoemission data to argue that quasi- 
particles are very strongly scattered over most of the 
Fermi surface. The only long-lived quasiparticles are 
found at the zone diagonals, where they decay with a 
weak T 2 Fermi liquid form. This is known as the 'cold 
spot' model and has been used to understand a wide Vja=. 
riety of experiments Csee for example Van der Marelta 
and Xiang and HardyEj). 

In this paper we consider this third scenario of 
anisotropic scattering and treat a minimal model which 
we believe is illustrative of the strengths and weaknesses 
of this approach (similar in spirit to, some treatments of 
the cuprate superconducting stataH3). Whilst more so- 
phisticated treatments exist, either with a specific re- 
laxation mechanismlljO or in— the limit of extremely 
anisotropic in-plane scattering!!!, there are advantages 
to this model's simplicity. We can obtain analytic ex- 
pressions for a wide range of transport properties over 
the full range of anisotropy. Our calculations illustrate 
the conclusions of Ioffe and Millis in the limit of strong 
anisotropy, but also allow us to see how this evolves from 
the more familiar isotropic metal. This is important as 
the approach to the isotropic limit may be more rele- 
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vant to experiments on overdoped materials where the 
two-lifetime behavior is less apparent. We also consider 
the transition to the intermediate and high field regimes 
of magnetoresistance within this picture, and non-linear 
effects at high electric field. 

Perhaps the key discrepancy between the 'cold-spot' 
model and the transport properties of the cuprates is the 
in-plane orbital magnetoresistance. The model predicts 
a large magnetoresistance with a distinct temperature 
dependence in contrast with current experiments— This 
was was noted in Ioffe and Millis' original papeitLj and 
is reproduced in passing here. It has been argued that 
proper inclusioruaf— vertex corrections could account for 
this discrepancycSEl The most important new work in 
this paper is to consider the orbital magnetoresistance for 
currents moving along the c-axis. The c-axis conductiv- 
ity in a quasi-two dimensional metal may be computed 
without vertex correctionso and so should be a robust 
feature of this scenario. We show how an in-plane mag- 
netic field affects the c-axis conductivity and argue that 
this experimental geometry provides a key 4-est of the 
model. Indeed experiments exist on Tl-2201t3 but only 
in the overdoped regime where the two lifetimes become 
less distinct. Nevertheless our complete solution of an 
anisotropic scattering model allows us to obtain a reason- 
able fit to the data provided we include both scattering 
anisotropy and bandstructure effects. 

Our main conclusions are that within our model the 
in-plane magnetoresistance remains an outstanding prob- 
lem if this model is to fit the experiment — not just from 
its magnitude and temperature dependence but from the 
scale at which deviations from the weak-field regime oc- 
cur. While a proper treatment of vertex corrections may 
correct this for in-plane properties, the c-axis magnetore- 
sistance is insensitive to these corrections and provides a 
robust test of the model. Experiments here have mainly 
been done on overdoped cuprates and our detailed fit 
to experiments are consistent with this model. High 
field c-axis measurements on optimally doped cuprates 
should be made to test unambiguously this transport 
phenomenology. 

This paper is planned as follows. In section II we 
present our models of scattering and bandstructure which 
are the ingredients for the analytic calculation of trans- 
port properties. The results are to be found in Sec- 



tion [II. These are displayed in the general case, and 
also in the limit of large scattering anisotropy. Both in- 
plane and out-of-plane magnetoresistivity are shown, and 
compared with experimental observations in section IV. 
Non-ohmic in-plane conductivity is also calculated. Con- 



clusions are drawn in section ^ 



II. ANISOTROPIC TRANSPORT 

Our main assumption in this work is to use a Boltz- 
mann equation within the relaxation time approximation 
as a description of transport. We write the collision in- 



tegral, /[/], as — rk[/(k) — /(°)(k)] and impose a phe- 
nomenological scattering rate, T^, which varies smoothly 
around the in-plane Fermi surface and has square sym- 
metry. 

We do not speculate on the underlying cause of the 
quasiparticle relaxation rate but, if this scattering were 
due to scattering from a soft bosonic mode at finite Q 
(antiferromagnetic spin fluctuations for example), then 
the validity of tLe relaxation time approximation might 
be questionableES. This is because such scattering would 
tend to equilibrate quasiparticles only between points 
on the Fermi surface connected by the spanning vec- 
tor (Q) of the mode. By contrast, the relaxation time 
approximation assumes that the quasiparticle is equili- 
brated uniformly around the Fermi surface. This will 
not affect our calculation of c-axis response if the bosonic 
mode is 2-dimensional, but it could change the in-plane 
physics. This is presumably the reason why vertex cor- 
rections could be important for in-plane transport. How- 
ever, including this process correctly would tend to make 
the quasiparticle distribution even more anisotropic when 
driven out of equilibrium. As we will see, the anisotropy 
of the out-of-equilibrium distribution of quasiparticles is 
precisely what is probed by the in-plane magnetoresis- 
tance. Increasing this anisotropy is likely to increase fur- 
ther the in-plane magnetoresistance. 

Within the relaxation time approximation we now 
introduce two forms of scattering rate to model the 
anisotropy. 



A. The cold-spot model 

Following Ioffe and MillisEI, the angle-resolved pho- 
toemission results suggest that there are points on the 
Fermi surface where the electron is a fairly long-lived ex- 
citation: points where the Fermi surface intersects the 
Brillouin zone diagonal. Elsewhere on the Fermi surface, 
only broad features in energy are seen in the photoemis- 
sion. Within a Fermi liquid picture one would then infer 
that electron-like quasiparticles are strongly scattered ev- 
erywhere except for these special points. We can capture 
this by the following form for our scattering rate 



r(fc||,6»,fcjj = T/cos 2 26» + T s sin 2 29 . 



(1) 



Here 9 is the azimuthal angle between the in-plane mo- 
mentum wavevector fcii and the a-axis direction of the 
crystal, as shown in Fig. |]. This is a simple model 
which nevertheless allows us to explore the transition 
from isotropic (r^ = T s ) to extremely anisotropic scat- 
tering (Tf r s ). The latter limit corresponds to the 
case studied by Ioffe and Millis, where, at the planar 
zone diagonals [9 = (2n + the scattering rate has 

a minimum. In their model, T s is just the Fermi-liquid 
scattering rate, T s = 1/tfl ~ T 2 . The scattering around 
the Fermi-surface in the limit of strong anisotropy is il- 
lustrated in Fig. |2|. 
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(a) 



(b) 



FIG. 1: The cylindrical coordinates used in this paper. This 
also illustrates the Fermi surface we consider with the c-axis 
dispersion greatly exaggerated. The degree to which the c- 
axis dispersion depends on the in-plane momentum varies 
within this model from (a) no dependence (7 = 0) to (b) 
the case when there is no dispersion along the zone diagonals 
( 7 = 1). 



In this form a — gives us the isotropic scattering case 
whilst a = 1 or — 1 gives large anisotropy with cold spots 
on the zone diagonals or zone axes respectively. We will 
express results in cither notation (Ty,r s or a, To) de- 
pending on applicability to available data and ease of 
interpretation, and are always able to calculate across 
the full range of a. 



B. The hot-spot model 

Alternatively it is has been proposed that hot spots 
exist near the (tt, 0) points of the Brillouin zone where 
strong scattering occurs, perhaps due to antiferromag- 
netic spin fluctuations. This leaves quasiparticles on the 
rest of the Fermi surface with a longer lifetime. To cap- 
ture this within our phenomenology we add the scattering 
lifetimes around the Fermi surface so 



r(Jfe|| , 9, k x ) = r f cos 2 26» + t s sin 2 29 . (5) 



Here we expect rt to be smaller than t s so that we have 
a model with strong-scattering hot spots and extended 
cold regions with less scattering (see Fig. ||) . 
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FIG. 2: The angular dependence of the scattering rate in 
our parameterization of the cold-spot model (dotted line and 
Eq. 0) and the hot-spot model (solid line and Eq. ^]). Although 
the current is preferentially carried by the longest-lived quasi- 
particle, in the cold-spot model the short-lived particles can 
contribute since they cover most of the Fermi line. In the 
hot-spot model the cold regions short circuit all transport 
processes and dominate all quantities. 



It will often be convenient to express the scattering 
rate in terms of an anisotropy parameter, a, where Eq. |l| 
is written as 



with 



T(k h e,k ± ) = r (l + a cos 40) , 



r = (r / + r s )/2 
.. _ (T/-r.) 



(T/ + r.) ■ 



(2) 

(3) 
(4) 



C. Bandstructure 

The next component of the model is the bandstruc- 
ture. We assume that the in-plane dispersion is isotropic, 
which is not unreasonable for the cuprates. By con- 
trast the dispersion in the out-of-plane direction is known 
to h | ay . o .s ig nificant dependence on the in-plane momen- 
tumE3E§E3. This is due to the local chemistry of the 
copper-oxide planes as emphasized by XiangEEI. Trans- 
port between planes occurs principally through the cop- 
per 4s orbitals. Since these orbitals are circularly sym- 
metric in the plane, their overlap amplitude with the 
Cu 3d x 2_ y 2-0 2p hybrids has d-wave symmetry as illus- 
trated in Fig. |^. The resulting tunneling term is therefore 
proportional to the square of this amplitude and can be 
represented by the form t± cos 2 (29). The most impor- 
tant effect this has on c-axis transport stems from the 
fact that the longest-lived quasiparticles on the zone di- 
agonals are precisely those that have a vanishing proba- 
bility of c-axis hopping. For this reason we ignore other 
features of the c-axis dispersion which may lead to ad- 
ditional small tunneling probabilities at other points on 
the Fermi-surface such as those due to a 
tetragonal unit cell considered elsewhereE 
structure is then given by 



dy-centered 
'. Our band- 



e(k||,0,k|_) =e(fc 2 )-2<_ L (l+ 7 cos46<)cos(fc i c). (6) 

To illustrate the importance of this effect and to allow 
for other tunneling mechanisms we parameterize the de- 
gree of tunneling anisotropy by 7, which we allow to vary 
from to 1. When 7=1, we have complete elimination 
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Cu 4s orbital 



FIG. 3: The overlap between the hybrids of copper 3d x 2_ y 2 
and oxygen 2p orbitals and a copper 4s has d-wave symmetry. 
Since transport between layers occurs primarily via the copper 
4s orbitals, the c-axis tunneling probability vanishes when the 
in-plane quasiparticle momentum is directed along the 
diagonals — the zeros of the d x 2_ y 2 form (after Xiang et al 



of c-axis transport for quasiparticles with in-plane mo- 
mentum along the zone diagonals. The frequency depen- 
dence of this end-point of the model has been studied by 
van der MareEj. Allowing both the degree of scattering 
anisotropy (a) and c-axis dispersion (7) to vary across 
their full range turns out to be crucial in understanding 
the out-of-plane transport and is a novel aspect, of this 
work. It unifies various other approachesE3EZrE3 which 
focus on the extremes of one or other of the ranges of a 
and 7. 



III. RESULTS 

Using the standard representation of the quasiparticle 
distribution, / = /(°) + ^<9 £ / (0) , we write the Boltzmann 
equation (for uniform fields and temperature gradients) 
afi 



dtij) + r k ^ - -r (E + v k x B) • V k ^ = 
n 

eE-v k - 



T 



v k -M7) 



We will solve this equation in various limits in the subse- 
quent sections. We will only be interested in d.c. prop- 
erties in this paper, however our results for the cold-spot 
model (Eq. [l]) are easily extended into the frequency do- 
main. After Fourier transforming with respect to time, 
we see that 



: 26» + r. 



'■29), 



-> (T f + iuj) sin 2 26 + (T s +iuj) cos 2 29 . 

(8) 

We can therefore obtain the finite frequency results for 
the cold-spot model from our expressions at d.c. simply 
by the prescription 



r / -> r / + iuj 



r. 



no 



(9) 



A. In-plane transport in weak magnetic fields 

We first focus on leading order response in an in-plane 
electric field. For a circular Fermi surface we may write 
Eq. as 

T(6)iJj - lu c ^ = eEv F cos6 , (10) 

where uj c — eBvp /(HIcf) and vf is the in-plane Fermi 
velocity. 

We may solve this equation by Jones-Zener expansion 
for in-plane currents with a weak magnetic field along the 
c-axis. The currents, and hence conductivity elements 
are calculated using the relation 



J/' 



eE 



V F 



ipdOdkj 



(11) 



where current is flowing in direction \i with velocity v^. 
We expand a xx and a xy to order B 2 and the resistivity 
is then found by simple matrix inversion. 

If we consider first hot-spot scattering (Eq. ||), we can 
see that this scattering rate can be eliminated as a viable 
model for the cuprates. Each of the properties evaluated 
is dominated by the long scattering time, which is the 
short-circuiting effect proposed by Hlubina and Ricetj. 
The in-plane conductivities are given by 



e 2 k F VF{Tf + t s ) 
Anhc 



(12) 



<7 (0) 
v xx 



u X X 



(3t 2 + 2t s t s + 3t 2 ) 
tan H = 0Jc f — ; — — + 0[B 3 } , 



4(T/ + T S ) 



, (2lT? - MT f T s + 2lT 2 



(13) 

0[B 4 } . (14) 



Within this model, all computed quantities are domi- 
nated by t s in the limit of strong anisotropy (r s 3> Tf). 
This model is unable to reproduce the two-lifetime phe- 
nomenology of the cuprates so we will not consider it 
further. 

However, for the cold spot model we obtain the follow- 
ing in terms of F f and T s : 



r(°) - 



e 2 vpk} 



1 



2^c ^/TfTs 



u xy 
yJ xx 



= tan Qh = — I =; — h — 



1 



1 



/ 



r. 



0[B 3 



(15) 
(16) 
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Cold-spot model Experiment Reference 
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TABLE I: A comparison of our results with experiment in the 
limit of highly-anisotropic scattering. 



an even more dramatic effect on the magnetoresistance 
which is higher order in magnetic field. We do not at- 
tempt to offer a solution for this puzzle though given the 
systematics in transport behavior across a wide range of 
cuprates it seems unlikely that this is due to fine tuning.. 
(By contrast in the charge- conjugation phenomenologyu 
the magnetoresistance is proportional to the square of 
the Hall angle and is small since this model is isotropic 
around the Fermi surface.) 
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0[B 4 ) . 
(17) 



This results in an in-plane orbital weak-field magnetore- 



sistance 



&Px 



o (0) 

\ J XX 



, .2 ( r / -r s ) 2 , 2 

Uc 81*1* (5r / 



7i>r s + 5r^) 



(18) 



which vanishes in the isotropic limit (r s — Tf). This 
reflects the well known result that there is no orbital 
magnetoresistance in an isotropic metal. In the high 
anisotropy (a — ► 1) limit, where Tf 3> T s , we find the 
results shown in table |. 

These results illustrate the way in which the cold-spot 
model reproduces the phenomenology of the cuprates as 
shown by Ioffe and MillisO. They take Tf to be large and 
temperature independent, while T s is assumed small and 
proportional to T 2 . The geometric mean then gives the 
linear in temperature resistivity while the inverse Hall 
angle is proportional to T s and hence varies as T 2 . The 
results are also illustrative of the problems of this model. 
Adding Zn impurities to the CuC>2 has been interpreted 
as adding a unitary scatterer which should add a temper- 
ature independent term to the scattering rates. Adding a 
constant elastic scattering term to the two rates appear- 
ing in the geometric mean which forms the resistivity 
(Eq. |l|) will not reproduce the Matthiessen's rule be- 
havior seen in experiment]]. Instead it will change both 
the temperature dependence and the intercept of the re- 
sistivity. 

The second difficulty is the large magnetoresistance 
that this model predicts. While experiments suggest that 
the magnetoresistance is proportional to the Hall angle 
squaredQ with a constant of proportionality which can 
be of order 1, the cold-spot model predicts that these 
two quantities differ by Tf/T s — a large, temperature de- 
pendent factor. It is tempting to appeal to new physi- 
cal processes to explain this (such as vertex corrections) . 
However the cause of this large magnetoresistance is in- 
timately related to the appearance of distinct tempera- 
ture dependences of the resistivity and inverse Hall angle. 
Magnetic fields couple to the anisotropy of ip around the 
Fermi surface as can be seen from Eq. [U] and so at each 
order in the magnetic field, the effect of the anisotropy 
in the model is amplified. This allows the inverse Hall 
angle to differ from the resistivity but consequently has 



B. Beyond the weak- field regime 

Thus far, in keeping with other workers, we have shown 
that in-plane transport is governed by a different scatter- 
ing rate (formed from the combination of Tf and r s ) at 
each order in the magnetic field. Here we go beyond the 
weak-field regime and show numerically that no new com- 
binations emerge beyond second order in the magnetic 
field until the magnetoresistance saturates. The scat- 
tering rate combination for weak-field magnetoresistance 
provides the scale for magnetic field effects throughout 
the intermediate field region. 

The saturation field may be obtained from a high field 
expansion of the Boltzmann equation — expanding the so- 
lution in powers of l/ui c . We find that the magnetoresis- 
tance saturates at a value of 



A Pxx (B -> oo) = (y/rj- VTT) 1 
Pxx(B = 0) 



(19) 



The high field limit is reached when w c ~ Ty so the 
fastest rate controls the transition to the high field limit. 
This is because saturation requires that quasiparticles 
can completely orbit the Fermi JMJrface and it is the hot 
regions which limit this processES. This is also the scale 
for Landau quantization. The approach to the saturated 
value is shown in Fig. ^. The magnetic field is plotted 
as a function of lu c /Tq since in the cold spot model this 
is dominated by Tf and so shows that this sets the scale 
for saturation. 

[As an aside, we can obtain the universal form of the 
magnetoresistance in an isotropic metal. Of course, this 
is formally zero but by dividing by the saturated value in 
the limit of small anisotropy we can obtain a finite result. 
This is also shown in Fig. [|.] 

While we cannot solve Eq. [l0| analytically except in the 
low and high field limits, the equation is readily solved 
numerically. The fact that we require a solution of ip 
which is periodic around the Fermi surface means that 
we can look for a Fourier series solution. Our numerical 
results reproduce the analytical results of the previous 
section. In addition, we find the remarkable result that 
the magnetoresistance is a function of a single dimen- 
sionless parameter until it saturates. This scaling pa- 
rameter governs the magnetoresistance even beyond the 
weak-field (B 2 ) regime of the magnetoresistance so we 
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FIG. 4: The in-plane magnetoresistance is shown beyond the 
weak field region as determined from a numerical solution. 
Saturation occurs at a scale determined by lo c ~ Fn since 
quasiparticles must precess around the whole Fermi surface. 
The inset shows that no new field scale emerges between the 
weak-field and the saturated regime. By scaling the data to 
the weak- field region (thereby defining r mr ) we see that the 
magnetoresistance continues to fall on a universal curve be- 
yond the B 2 regime until the saturated magnetoresistance 
value is reached. 



may write 



Ap xx (B) 
Pxx{0) 



lu c < r 



(20) 



where r mr is defined by the weak-field magnetoresistance 
such that Ap xx /p xx = (w c r mr ) 2 . For Tf ^ T s we have 

tL ~ ^ f /8T 3 . 

We illustrate this universality in the inset to Fig. 0. 
We see that as a function of this combination of scatter- 
ing rates, the magnetoresistance follows a universal curve 
until the saturation value is approached. An intermedi- 
ate field region emerges for fields 1 < cj c T mr < w c /T/, 
where the magnetoresistance adopts a new power law 
~ B 33 . Thus our numerical treatment suggests that 
no new scattering rate combinations are produced in the 
Jones-Zener expansion until the magnetoresistance sat- 
urates. We have confirmed this is the case analytically 
at fourth order in the magnetic field. This scaling also 
means that in our model there is an absolute value of 
magnetoresistance above which one must see a deviation 
from a B 2 dependence. Our model predicts that with 
a magnetoresistance of only 0.5% one would expect at 
least a 10% deviation fmm a B 2 form. Comparing this 
with experimental workH on optimally doped Tl-2201 in 
pulsed magnets, we find a deviation from quadratic de- 
pendence is only perhaps being seen in optimally doped 
Tl-2201 at a magnetoresistances of order 3%. (The ob- 
served deviation from the weak-field region is more con- 
sistent with the Coleman-Schofield-Tsvelik phenomenol- 
ogy which predicts that this should occur when the mag- 



netoresistance is at or less than about 109cEj.) 



C. c-axis transport 

In-plane magnetoresistance in the cold-spot model is at 
odds with the current experiments. It has been ar. 
that including vertex corrections could correct thisESI 
As we have argued in the introduction, c-axis transport 
is then of fundamental importance in testing the cold- 
spot model. This is because to lowest order in t\ the 
c-axis conductivity is related to the convolution of two 
in-plane spectral functionsEJ. These are the quantities 
which may be inferred from the angle-resolved photoe- 
mission that inspired the cold-spot model. Thus, inde- 
pendent of a picture of quasiparticles for in-plane trans- 
port, a cold-spot model should be able to account for the 
c-axis d.c. conductivityE-1 In particular we must show 
that this leads to a consistent picture of c-axis orbital 
magnetoresistance (magnetic field in the plane and elec- 
tric field out of the plane). The most systematic Sjfcqdy 



of these effects has been performed by Hussey et altB in 
overdoped Tl-2201 and this motivates the following anal- 
ysis. The two-lifetime behavior is much less apparent in 
the overdoped material so we will not be able to use the 
Tf 3> T s asymptotics. Instead we need a full solution for 
any degree of anisotropy. 

We will now calculate the c-axis magnetoresistance, 
Ap zz / p zz . To obtain this, we first calculate the relevant 
components of the conductivity matrix, Uij and invert. 
In this experimental geometry, a has zero components 
Cxy,cr XZl iJy X and a zx due to the B field being in-plane. 
a xx is equal to the zero-magnetic field value of a yy . We 
can expand the other terms in magnetic field: 



71 — 3 

_ (0) , V- (n) 



(21) 



where the superscripts refer to the order of effect in mag- 
netic field. Symmetry under time reversal means, of 



course, that a zz 



and a z 3 J 



and the Hall terms a y °J and 



T (2) 

> yz 



are zero. 



In addition the Hall term Oy 1 } is small 



(~ t\). Under these conditions, the c-axis magnetoresis- 
tivity in weak magnetic field simplifies to 



0[B 6 



Pzz 


= pi°J + 




Api 2 } 


a?} 




pi } 


V zz 






( (2 
(Viz 


V zz (J zz J 






CT (0)2 

V zz 



(22) 



He*e we have adopted tha-sign convention of Hussey et 
alxB and Dragulescu et alEi, expecting the fourth-order 
magnetoresistivity to be negative. 

To obtain the conductivity, we follow the same general 
procedure as for the in-plane case, solving the Boltzmann 
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equation using a Jones-Zener expansion. We now intro- 
duce the angle, cj>, which is the in-plane angle between 
the B field and the a-direction. becomes the azimuthal 
coordinate relative to the direction of the B field. The 
solutions to the Boltzmann equation are now 



eE 



dc 



hT(0 + (j>) dk ± 
where, to lowest order in t±, 
eB n (\ 



X> (n) (k), (23) 



^ (n) ( k ) 



TiT{6 



-v F sinf— 

n ok± 



(24) 



The ratio of azJ/azJ, when expanded in powers, of a 



yields a zero-order term which is, as Hussey et alMB first 
showed 



a (2) 

f zz 
Ozz 



•Of 



(25) 



Similarly, for the quadratic term in the c-axis conductiv- 
ity: 



cjj ( a K zV ^ _ 6 + 3 7 2 + 2 7 cos(4 </>) 
(0)'\J0)) - 2(2 + 7 2 ) ' [ ' 



.(2). 



a 



•a. 



which in the limit of simple inter-plane tunneling (7 — * 0) 
gives 3/2, again as shown by Hussey et al. 

The Hall conductivity, a yz , is found to vary as t\. In 
the limit of isotropic in-plane scattering the out-of-plane 
Hall angle is identical to the in-plane oneo 



lim a _ 



cBvf 



A ) hk F T 



(27) 



This is not true in general however. Nevertheless its de- 
pendence on tj_ means that the Hall conductivity does 
not appear, at leading order, in the out-of-plane magne- 
toresistance. 

Rather than display the full expression for every quan- 
tity we have calculated, we show graphically two quan- 
tities which Hussey et al. examined experimentally in 
single-crystal Tl2Ba2CuC>6. There a key observation was 
the four-fold variation of the c-axis magnetoresistance as 
the magnetic field was rotated in-plane, i.e. that 



(28) 



where p~tz and ffz} were both found to be positive. In 
fact it is straightforward to show that this 4-fold modula- 
tion at fourth order in B is a simple consequence of square 
symmetry. Hussey et al. initially analyzed their angle- 
dependent magnetoresistance results purely in terms of 
anisotropic, in-plane scattering (as parameterized here 
by a). We first look at the offset part in the fourth order 
magnetic field term in the magnetoresistance. We exam- 
s' 4 ' /o' ' 

- — a quantity found experimentally to be 
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and (b) experimental datac3 from 



measurements on Tl2Ba2Cu06, where this quantity was 0.6± 
0.1, approximately independent of temperature. This implies 
that a should not vary greatly in the experimental range of 
temperatures. 
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FIG. 6: (a) The size of the fourth order modulation in c- 
axis magnetoresistivitpTirelative to the offset at this order and 
(b) experimental dataEj from measurements on Tl2Ba2CuOe 
again shown for comparison. A positive value of fizz lpzz 
means that resistance maxima exist for a magnetic field along 
[110]. The effect of bandstructure is therefore pronounced and 
produces resistance maxima where the cold spots are on the 
zone diagonals (a > 0). In (a), the approximate constraints 
on the valid regions of (a, 7) parameter space are given by the 
limits of the modulation as taken from experiment. These are 
shown as dotted lines at values of 0.15 and 0.05. 



roughly independent of temperature and approximately 
equal to 0.6 (0.6 ± 0.1). We show it for various values of 
7, as a function of a in Figure [|. This puts a constraint 
on the temperature dependence of a, an issue which is 
explored in section [V where we try to fit magnetoresis- 



tance data to our model. 

The second quantity studied experimentally is the am- 
plitude of the modulation in magnetoresistance at order 

B 4 . The quantity evaluated is pit /pit , which is pos- 
itive when the resistance is maximum with B aligned 
along the zone diagonal (Eq. ^8|). This too is plotted 
across the range of scattering anisotropy, a, for various 
bandstructures, 7, in Figure ^. 

Hussey et al. found it to vary with temperature, drop- 
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ping from about 0.15 at 50K to 0.05 at about 140K, and 
these limits are shown as dotted lines in Figure [j. In 

our calculations, pzz/pzz is positive when a is nega- 
tive, corresponding to a situation where the hot spots 
and cold spots have effectively exchanged their usual po- 
sitions. This can be explained intuitively by identifying 
the maxima in the Lorentz force, v x B, at points on 
the Fermi surface orthogonal to the in-plane magnetic 
field. Then, loosely, the parts of the Fermi surface least 
affected by a magnetic field are those with Vp parallel to 
B. Thus, naively, one would expect the minima in the 
magnetoresistance to occur when the field is along the 
cold-spot directions. In fact experimentally this orienta- 
tion gives a maxima in the magnetoresistance. Turning 
on anisotropy in the c-axis dispersion 7, one can suffi- 
ciently inhibit c-axis transport through the cold spots in 
the positive a regime that resistance maxima can be seen 
for B along [110] with cold spots along [110] in keeping 
with the experiments. _ . 

This effect has been attributed elsewhereEil to the cur- 
vature of the real Fermi surface in a model of isotropic 
scattering. Here we have shown that the two effects will 
work together in the regime of positive bandstructure 
anisotropy. Once again, however, we see that there is a 
constraint on the allowed degree of variation of scatter- 
ing anisotropy with temperature. In Table [n] we show 
the form of these 4th order quantities shown in Figures 
H and H as well as the size of the 2nd order magnetoresis- 
tivity and in-plane properties. This will be useful when 
fitting to data later. 



which all simplify or reduce to zero appropriately in the 
isotropic limit, a — > 0. Furthermore, in the limit of 
strong anisotropy, jO,— > 1, we find o\ — — 02, as demon- 
strated by Hlubinatj. 

The most obvious consequence of non-linear response 
in the electric field is that the current no longer flows 
parallel to the electric field. Defining the parallel current 
(j||) as being the response parallel to the applied field, 
and the transverse current as the current component per- 
pendicular to the applied field (J±) we find that 



3-L 



-E 3 (a 1 — a 2) sin^ 



(35) 



where <p is the angle between the in-plane electric field 
and the a-axis. Using our expressions for o~\ and o~2 
(Eqs. p3 and B4T) we may write 



3± 

3\\ 



( eE 



lim 



a (2 + a 3 ) 
2(1 -a 2 ) 3 



(36) 
(37) 



Here Vd = eE/2mTo is the Drude drift velocity of a fast 
decaying quasiparticle. So, although the Vd/vF is gener- 
ally tiny, nonlinear effects are indeed strongly enhanced 
by anisotropy. 



E. Thermal transport 



D. Non-ohmic conductivity 

We have also calculated the non-ohmic, in-plane coh-i 
ductivity of the system in zero magnetic field. HlubinaE£l 
has suggested that this is a useful test of the model of 
Ioffe and Millis. The Boltzmann equation (Eq. ^) is now 
solved by Jones-Zener expansion in E and is given by 



TkV' - r E • 



VkV' = eE • v k . 



(29) 



From the 4-fold symmetry of the system, it is straightfor- 
ward to deduce that the currents in the x and y directions 
must have the following form at third order in the electric 
field: 



Jx 

jy 

We then find that 



a E x + aiEl + a 2 E 2 E x , 
a E y + a x El + (r 2 E 2 E y . 



tfi = - 



e 2 kpvp 



0~2 = 



2TThcTo(l -a 2 ) 1 / 2 ' 

e 4 v F a((l + a)(a 2 +2a + 2)) 
^Ti 3 k F cT 3 (l-a 2 ) 7 /2 
e A v F a((a 3 - 9a 2 + 4a - 6) 



(1-a 2 ) 7 / 2 



(30) 
(31) 

(32) 
(33) 
(34) 



In this treatment we have not speculated on the 
anisotropy of the energy relaxation rate around the Fermi 
surface. However if we make the assumption that the en- 
ergy relaxation follows the quasiparticle relaxation rate 
and only varies around the Fermi surface (as opposed to 
variations away from the Fermi surface) we can calculate 
all of the thermal transport coefficients. Without repro- 
ducing the details of the calculation we can see immedi- 
ately from the form of the Boltzmann equation (Eq. [7]) 
that temperature gradients drive the quasiparticle distri- 
bution in exactly the same way as electric fields except 
for the usual factor of (ek — /J,p)/T. This is seen in the 
right-hand-side of Eq. 0. This factor is isotropic around 
the Fermi surface and does not introduce any further an- 
gular dependence. Thus we can conclude that all currents 
proportional to temperature gradients will have exactly 
the same dependence on scattering rate as those driven 
by electric fields (for example yT fT s in the absence of a 
magnetic field). Within this approximation the thermal 
conductivity will obey the Wiedemann- Franz law and the 
diffusion thermopower will have the usual linear temper- 
ature dependence and be independent of the scattering 
rate. To account for the unusual systemalics in the mea- 
sured thermopower of the cuprate metala^J in this model 
one would need to add a significant anisotropic energy 
dependence of the scattering rate that differs from the 
anisotropy of the transport relaxation rate. 
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TABLE II: Results of our calculations of in-plane and c-axis properties. In particular we have attempted to find dimensionless 
quantities, fx, fi and fz, which are functions of a and 7. This enables us to constrain the values of a and 7 by comparing 
directly to experiments without further assumptions. 



IV. APPLICATION TO EXPERIMENTAL DATA 

In order to apply the analytic results of our model to 
experimental data, we have made a comparison of mag- 
netoresistivity results with measurements on overdoped. 
samples of single-crystal Tl 2 Ba2Cu06, made by TylerEa 
and Husseyc3. We choose this system because its c-axis 
magnetoresistance is well characterized. We have argued 
that c-axis properties are the most robust quantities in 
this model since they are not affected by vertex correc- 
tions, nor do they rely on a quasiparticle picture. We 
then identify combinations of measured c-axis quantities 
which directly probe the degree of anisotropy with min- 
imal dependence on unknown parameters. Finally we 
address the in-plane transport features. We also iden- 
tify combinations of experimental quantities which allow 
a comparison with theory involving the fewest assump- 
tions on unknown parameters. In particular we can con- 
strain the degree of anisotropy in the c-axis dispersion 
(7) and the range of scattering anisotropy (a) from the 
c-axis magnetoresistance. 

First we consider the overall magnitude of the fourth 
order magnetoresistance by comparing this to the sec- 
ond order term. This is shown in Fig. |E|(b). When we 
compare this to our model [Fig. ||(a)] we see that to ob- 
tain the observation of an essentially temperature inde- 
pendent result of around 0.6, the degree of scattering 
anisotropy must be limited to a < 0.5. Furthermore, if 
a is to have any temperature dependence at all then 7 
should be approaching 1 where the gradients in Fig. |^(a) 
are smallest. 

We have already remarked on a second conclusion from 
a comparison between theory and experiment. In order 
for the maxima in the c-axis magnetoresistance to occur 
when the field is along the zone-diagonals, c-axis trans- 
port along these directions must be suppressed. This is 
illustrated in Fig. ^(a) where we see that with a > we 

require 7 > to obtain a ratio of p^J /pzz of the correct 
sign. Furthermore the experimental bounds on this ratio 
[illustrated in Fig. ||(b)] confirm that the degree of scat- 
tering anisotropy (a) must be less than about 0.6. In no 
sense then are we in the limit of strong anisotropy in this 
overdoped material. 

A further observation comes from the temperature de- 
pendence shown in Fig. @(b). Generically the quantities 
shown in Figs. |5|(b) and g(b) should be temperature de- 
pendent so the observation that only the second of these 



has any significant variation with temperature restricts a 
to vary between about 0.2 and 0.4 with a band structure 
anisotropy 7 ~ 0.9. Also, in order to fit the trends in the 
data, we see that the overdoped state must become less 
anisotropic as the temperature is lowered. Proximity to 
a zero temperature critical point would suggest the op- 
posite should be true but perhaps here elastic processes 
are begining to dominate. 

In order to predict other experimental quantities, we 
then need to know the temperature dependence of To. 
This too can be found from the second order c-axis mag- 
netoresistance, when data is compared with the func- 
tional form in table |n|, here fixing 7=1. In addition 
it can also be independently infered from the in-plane re- 
sistivity but now relying on the assumptions of Boltzman 
transport theory. 

With these two temperature dependencies known, we 
may try to predict in-plane properties. We show plots of 
the cotangent of the Hall angle, cot 0#, and the in-plane 
magnetoresistance ^.pxx/pxx, together with experimen- 
tal data from Hussey and Tyler in Figs. ^ and 0. We 
see by comparing the experimental data with the param- 
eters extracted from other measurements that the data 
and theory follow the same trends and are of the cor- 
rect order of magnitudes. The two different fits on each 
plot represent the two different methods of infering Tq: 
either using the in-plane resistivity or the out-of-plane 
magnetoresistance. The difference in these two methods 
is that the in-plane measurements would be expected to 
be sensitive to the transport lifetime, as opposed to the 
quasiparticle lifetime in the case of out-of-plane measure- 
ments. Surprisingly it appears that the quasiparticle life- 
time is longer than the transport lifetime. 

Overall we see that for overdoped Tl-2201 a consis- 
tent picture emerges of a Fermi-liquid-like metal with 
weak scattering anisotropy around the Fermi surface but 
a strong anisotropy in the c-axis dispersion. A more 
rigourous test of the model would be obtained by per- 
forming a similar comparison on optimally doped Tl-2201 
(T c = 85K). There the two-lifetime behavior of in-plane 
magnetotransport is clearly apparent with vsey good fits 
to data being found with simple power lawscS: 

p xx = -20 + 1.56(T/K) (iHcm , (38) 
(B/T)cot0 H = 300+ 1.78* 10~ 2 (T/K) 2 , (39) 

for 100K < T < 300K. In the limit of high in-plane 
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in-plane properties to predict the magnetoresistance so 
uncertainties about vertex corrections are present in this 
estimate. 

Again, using the in-plane properties, we would expect 
that the out-of-plane magnetoresistance would be 0.1 in 
fields as low as 3T with a fourth order magnetoresis- 
tance of 0.001 in an in-plane field of about 5T. These 
effects should be measurable but again rely on using in- 
plane transport lifetimes to provide a measure of in-plane 
quasiparticle lifetimes. However, we can make an unam- 
biguous prediction within this model: that the positions 
of the maxima in resistivity at fourth order should move 
from the [110] directions (seen in the overdoped system) 
to the [100] directions in the optimally doped materials. 
This is seen in Fig. ^|(a) where for large a we see that 
bandstructure is unable to prevent the sign of the mod- 
ulation term, pzJ /pi^z, from being negative. 



FIG. 7: A comparison of ourr- theoretical prediction for 
cot 8h with experimental dataEj from measurements on 
ThP^CuOe. Predictions using Fo found from both in-plane 
and out-of-plane properties are shown, relating to using the 
transport or the intrinsic quasiparticle lifetime. 
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FIG. 8: A comparison of our theoretical prediction for in- 
planp^i second order magnetoresistance with experimental 
dataa from measurements on Tl2Ba2CuOe. Again, predic- 
tions using To found from both in-plane and out-of-plane 
properties are shown. 



anisotropy we can combine these measurements to ob- 
tain the magnetoresistance. We find 



o {0) 

I 1 XX 



7T 2 C 2 



o {0) 

\ J XX 

cot 2 e H 



(40) 



This would predict a magnetoresistance of about 0.04 at 
10T and 130K (so already out of the weak-field regime) — 
a factor of 40 greater than currently seen. We have used 



V. CONCLUSION 

We have presented a thorough investigation of a min- 
imal model for transport in a quasi-2D system, where 
we allow for anisotropic scattering in-plane and an 
anisotropy in the out-of-plane dispersion, both fully vari- 
able. Such a model has been proposed for the normal 
state of the cuprate superconductors with either strong 
scattering hot spots and weak scattering cold regions 
around the Fermi surface or hot regions with cold spots. 
We have studied the transport properties of both types 
of models and illustrated how short-circuiting in the hot- 
spot model makes this inconsistent with transport mea- 
surements on the cuprates. For the cold-spot model we 
have computed in- and out-of-plane magnetoresistivity, 
in-plane non-ohmic conductivity and thermal conductiv- 
ities for an arbitrary degree of anisotropy. 

We find that the in-plane magnetoresistance in this 
model is too large when compared with experiments on 
the cuprates, in keeping with other workEH. In addition 
the in-plane magnetoresistance should be universal be- 
yond the weak field limit with a well defined deviation 
from a B 2 dependence at weaker fields than currently ob- 
served. It has been argued that vertex corrections may 
account for these discrepancies between the model and 
experiment. Here we have focussed on the out-of-plane 
magnetoresistance for which there are no vertex correc- 
tions in the quasi-2D limit. We have completely char- 
acterized the magnetoresistance to fourth order in this 
geometry. We have compared our model with experi- 
ments on overdoped T12201 and show that the model 
can be reasonably well fit to the experiments. This re- 
quires only weak scattering anisotropy but a high degree 
of bandstructure anisotropy in the c-axis dispersion. 

A better test of this model would be to compare it to 
optimally doped Tl-2201 where the two-lifetime behavior 
is very clearly seen in in-plane magnetotransport mea- 
surements. This model would predict that strong angle 
dependent c-axis magnetoresistance should be observed 
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and the positions of the minima in the optimally doped 
system should be rotated by 7r/4 from those of the over- 
doped material. A quantitative analysis of the degree of 
in-plane anisotropy could also be made. The extent to 
which vertex corrections control the in-plane transport 
properties could then be assessed. 
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